The typical behaviour of relays 
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The typical behaviour of the relay-without-delay channel and its many-units generalisation, 
termed the relay array, under LDPC coding, is studied using methods of statistical mechanics. 
A demodulate-and-forward strategy is analytically solved using the replica symmetric ansatz which 
is exact in the studied system at the Nishimori's temperature. In particular, the typical level of 
improvement in communication performance by relaying messages is shown in the case of small and 
large number of relay units. 
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I. INTRODUCTION 

Methods of statistical mechanics have recently become 
increasingly more important in the study of communi- 
cation channels. The development of the replica and 
cavity methods for analysing disordered systems 0, 0] 
and the related recent introduction of systematic rigor- 
ous bounds [H, 0| made new theoretical tools available for 
their analysis. 

More specifically, the replica method has been applied 
to a wide range of problems in information theory, from 
error correcting codes [H, [|| to multiuser communica- 
tion [7|. It facilitates the derivation of practical and 
theoretical limits in various communication channels and 
provides typical results in cases that are difficult to tackle 
via traditional methods of information theory. 

The growing use of information networks, both physi- 
cally connected and wireless, and the increasing number 
of services taking place in the Internet, have made the 
study of multiuser communication highly attractive and 
relevant from a practical point of view, in addition to 
being a challenging and exciting field for theoretical re- 
search. 

Up to date, there is no generalised theory of multiuser 
channels within the framework of information theory and 
analytical results are only known for special cases; the 
main difficulty being that multiuser networks do not ad- 
mit the source-channel separation principle which plays 
an essential role in the analysis of communication chan- 
nels. Nevertheless, multi-user communication plays an 
important role in a variety of communication devices 
ranging from mobile phones to computers. We strongly 
believe that a statistical physics-based analysis may offer 
answers where the current information theory method- 
ology fails, especially in the limit of a large number of 
users. 

With the technological demand and the possibility of 
a more thorough study by the methods of statistical me- 
chanics, early results for multi-user communication are 
being revisited and analysed from different and comple- 
mentary points of view, resulting in new insights and 
developments @, H|. One of those interesting types of 
channels is the relay channel @ . This generic channel is 



characterised by an auxiliary user between the transmit- 
ter and receiver, which helps in the transmission of the 
message. Due to the increase in the number of multi-user 
networks, like mobile phones and computer networks, the 
transfer of information with the help of a relays became 
an attractive option. As these networks are becoming 
more distributed, the transmission with the help of arrays 
of relays becomes feasible and merits further analytical 
exploration. 

This paper is organised as follows. In section |TT] we 
define the general relay array and introduce as particular 
cases the classical relay channel and the relay-without- 
delay. In section IIIII we outline the statistical physics 
methods used to analyse the problem which will be based 
on a replica approach detailed in section [IVJ Section [V] 
contains our conclusions and final comments. 



II. THE MODEL 
A. LDPC Codes 

Low-Density Parity-Check (LDPC) codes [l(| are 
state-of-the-art error-correcting codes with performance 
that is second to none, especially within the high code 
rate regime. In the notation we will be using here, N- 
dimensional messages s are encoded into Af-dimcnsional 
codewords t. LDPC codes are defined by a binary parity- 
check matrix A = [C\ | Ca], concatenating two very 
sparse matrices known to both sender and receiver: C2 
that is invertible and of dimensionality (M—N) x (M—N) 
and C\ of dimensionality (M — N) x N. The matrix A 
can be either random or regular, characterised by the 
number of non-zero elements per row (K) and column 
(C). Irregular codes show superior performance to reg- 
ular structures [ill, E3 if constructed carefully. In or- 
der to simplify our treatment, we focus here on regu- 
lar constructions; t he g eneralisation to irregular codes is 
straightforward [13.[l4|. 

Encoding refers to the linear mapping of a N- 
dimensional original message s S {0, 1}^ to a M- 
dimensional codeword t e {0, 1} M (M > N) 



t = Gs (mod 2) 



(1) 
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FIG. 1: The L-component relay array. The original sender 
sends a codeword t to the final receiver and to each of the L 
relays. Each relay receives a message r; which is a corrupted 
version of the original codeword subject to the AWGNs vo and 
Vi. It then sends to the final receiver the encoded vector ti. 
The final receiver receives the original transmitted codeword 
summed with all t,'s and corrupted by the AWGN vectors i/o 
and v. 



where all operations are performed in the field {0, 1} and 
are indicated by (mod 2). The generator matrix is 



G = 



(mod 2) 



(2) 



where I is the N x N identity matrix. By construction 
AG = (mod 2) and the first N bits of t correspond to 
the original message s. 

Decoding is carried out by estimating the most proba- 
ble transmitted vector from the received corrupted code- 
word [f| [l3|. For mathematical convenience, in the 
present work we map the Boolean variable t £ {0, 1} M 
into a spin variable t £ {1,— 1} M by the transformation 
x -> (-l) x . 



B. The Relay Array 

The relay array is a many-units generalisation of the 
(single unit) relay channel of 0. The LDPC codeword t 
is transmitted to each one of L relay units through noise 
channels corrupted by a global Additive White Gaussian 
Noise (AWGN) u and by local independent AWGNs i/,-. 
Each relay processes the received corrupted message 
and encodes the acquired information into a vector tj 
which is then transmitted to a final receiver. The final 
receiver receives an algebraic summation of the relay out- 
puts plus a direct transmission from the original sender, 
corrupted also by uq, subject to a final AWGN v. The 
exact form of the channel is depicted in Fig. [T] and the 
corresponding equations are 



r = at 



L 

£ 

i=l 



biU + v + v , 



at + v t + u . 



(3) 
(4) 



The variables a, bi and a (i = 1, L) are the relative 
gains of each transmission and can be random or set to 



constant values. The power from the original source to 
each relay is cj , to the final receiver is a 2 and the power 
from each relay to the final receiver is bf. 

When L = 1, we refer to the channel simply as the relay 
channel. In the classical relay channel (CRC), studied 
by Cover and El-Gamal [l5[, the messages sent by the 
relays to the final receiver are only allowed to depend 
on the set of symbols received by the each of the relays 
before the current time step, if = /(i 1 , i M_1 ), which 
corresponds to the fact that it takes the relay some time 
to process the information before relaying it. However, 
if the time delay in the direct transmission to the final 
receiver is much longer than in the transmission to the 
relay units, we can allow the message sent by the relay to 
depend on the present received symbol as well such that 
t% = /(i 1 , ...,#*). This last case, termed relay-without- 
delay (RWD), created significant interest recently and 
was studied by El-Gamal and Hassanpour [l6[. For the 
case of a relay array where all communication is carried 
out through the relays and there is no direct transmis- 
sion to the final receiver, the restriction of the CRC, to 
consider all but the last received symbol, is unnecessary. 



III. STATISTICAL PHYSICS OF DECODING 

We transform the decoding problem of the final re- 
ceiver into a statistical physics system by defining a dy- 
namical variable t = (t 1 ,...,t m ) which represents the 
candidate variable vectors at the receiver. Each plays 
a role equivalent to a spin located in the /i-th site of a 
lattice with M sites. 

The final receiver generates an estimate t of the orig- 
inal codeword using the Marginal Posterior Maximiscr 
(MPM) estimator 



< M = sgn (t 1 - 



(5) 



which minimises the probability of bit error [13|, 1 1 71 ] - 
Other estimators can be used depending on the error 
measure considered. For example, minimisation of block 
error is obtained using the Maximum a Posteriori (MAP) 
estimator. 

The posterior probability density is calculated by 
B ayes' rule as 



V(r\r) = 

with 

P(r) = X>Mt)P(t) 



nr\T)V(r) 
V(r) 



(6) 



t,{U},{r*} » 

(7) 

One of the basic quantities of interest is the overlap 
between the codeword and the decoded message. Our 



3 



analysis, focuses on the typical behaviour of the decod- 
ing process and, accordingly, we take averages over all 
possible codewords, all received messages and all allowed 
encodings, which we consider as quenched disorder in the 
system. The overlap between decoded and original mes- 
sage takes the form 



1 M 



(8) 



This quantity can be derived from the free-energy 



with the partition function 



and the corresponding Hamiltonian 

W(t;r) =-lnP(r|t)P(t). 



(10) 



(11) 



Usually we disregard the normalisation of the distri- 
butions within the Hamiltonian as they merely add con- 
stants that shift the zero energy. In the case of LDPC 
codes, 'P(t) turns out to be a constraint on the summa- 
tion variables. 

In the above Hamiltonian, the parity-check matrix A 
defines an interaction between the r variables while t 
and r act as local fields. The inverse temperature f3 is 
the ratio between the true and the decoder's assumed 
noise level. In our numerical calculations, we adopt (3 — 
1, also known as Nishimori's temperature, which means 
that the decoder assumes the correct noise level for the 
channel. It can be shown that at Nishimori's temperature 
the system never enters the glassy phase [2|, [l8| and the 
thermodynamically dominant solution is always Replica 
Symmetric (RS); we therefore restrict our analysis to the 
RS treatment. 

One of the important properties and the novelty of 
the statistical physics formulation of the problem is that 
looking at the problem as a dynamical spin system, we 
can interpret the results in terms of phase transitions 
which can be characterised by the overlap and the en- 
tropy function. Combining this extra information we 
can have a better understanding of the way the system 
changes from a phase of perfect decoding (which is called 
the ferromagnetic phase) to a phase where the message 
is recovered only up to a certain amount of error (the 
paramagnetic phase). 

The most studied strategies used by the relay units are 
the Amplify-and-Forward (A&F) and the Decode-and- 
Forward (D&F) strategies. In A&F, the relay just re- 
transmits its received vector, i.e., = Vi. As the replica 
treatment of this case turns out to be the same as for 
the simple Gaussian channel with a modified power and 
noise level, the solution is obtained straightforwardly by 
applying the results of [1 31 ] and will not be studied here. 



In the D&F strategy, the relays decode the message 
and transmit their estimate to the final receiver. Full use 
of LDPC decoding in the relays is made when each relay 
decodes the received vector by the MPM estimator us- 
ing the fact that the codeword was encoded by an LDPC 
code. The message transmitted to the final receiver by 
each relay would then be 



In equation ^ this is equivalent to setting 

M 



V(t i \r i )= JJj^-sgn. 



'i /•p(T«|r j 



^=1 



(12) 



(13) 



As t^ £ {±1}, we can rewrite this probability density 



as 



M 



(14) 



where 9 is the Heaviside step function. 

The replica treatment of the LDPC D&F turns out to 
be extremely involved due to the introduction of a theta 
function with an average over the variables inside it, 
which includes a term dependant on the parity-check ma- 
trix. Analytical studies in order to solve this, rather dif- 
ficult case, are under way. 

In the present work we focus on a simplification of this 
strategy also known in the literature as Demodulate-and- 
Forward. In it case, the relays do not have the complete 
information about the encoding mechanism and therefore 
assume a uniform prior for the transmitted codeword. In 
this case, the posterior distribution of the bits in the 
message for the relay is 



M 1 

'P(^i\ r i) = I I - : f „ u u I, 9 

1 + ex P h 2 vf 



(15) 



and it is straightforward to show that the MPM estimator 
is given simply by 



ti = sgn(rf). 



(16) 



The fact that the disorder relative to different codes 
does not appear in the estimate of the relays makes the 
replica calculations feasible in this follows. 



IV. REPLICA SYMMETRIC ANALYSIS 

As the RS analysis of LDPC coding systems has been 
introduced and carried out in a number of publications 
(e.g. [13j) we will skip the detailed derivation and con- 
centrate on the final expressions. The derivation follows 
exactly the same steps as in [l9| where quenched aver- 
ages are first carried out, followed by the RS assumption 
which enables the representation of the order parameters 
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in the form of field distributions. These are obtained us- 
ing a set of self-consistent saddle point equations of the 
form 



(17) 



= S x- 





£ T T[tt(r,r)] g n£ri 1 ( 1 + ™ 



(Tj - C£r) 

"2^ + ^) 



x exp 



Y2 bi s s nj 



(18) 



The expression (ti) is the mean of the variable U and 
V(r) oc ^(Ijt). The overlap is 



(sgn u) 



T{u) = { 5 [u- 



with 



£ r [*(r,r)^n^i(l + ^ 



(19) 



(20) 



the free energy is given by 

/3/ = -|ln2 + C(ln(l+xi)) 3 



*M ,+ fi-)). 

(ln|^[vI/(r,r)]' J n( 1 + ™0 



and the internal energy, the derivative with respect to /3 
of the above equation, is 



£ T ^(ln*)nti(l + ^ 
, Er^n^lCl+T*) 



(22) 



For any number L of relays, the results can be obtained 
by a numerical solution of the equations (fT7|) . Note the 
summation over the internal variables, i.e., the messages 
received and sent by the relays. This comes from the 
Bayesian formulation of the problem where the final re- 
ceiver has access just to r and, therefore, must integrate 
over the unknown variables. 

We also note that the above equations are fairly gen- 
eral. Using the appropriate "potential" \& we can recover 
all previous results for single user channels and apply 
them to more general channels when the intermediate 
processing of the message does not involve the previous 
knowledge of the parity-check matrices. 



The ferromagnetic state, which corresponds to perfect 
decoding, is given by the following solution to the saddle 
point equations (fTT|) 



(23) 



7r(i) = S(x — 1), 7r(x) = S(x — 1). 



Substitution of these distributions in equation (|T9|) 
gives d = 1. By substituting the ferromagnetic solution 
into the formulas for the free and internal energies, we 
obtain (for Nishimori's temperature) 



/ = -<ln#(l,r)) r , 



(24) 



meaning that the entropy of this phase is zero. 

The Hamiltonian of the relay array is gauge invariant 
with respect to the gauge transformation 



7^, 



(25) 



where the vector *y obeys the parity-check constraints. 
We can verify that the transition probabilities V{r^\t^) 
are also invariant under this gauge transformation. Note 
that if a channel is symmetric (for a definition see [I(|), it 
is automatically gauge invariant under the above trans- 
formation. For gauge invariant channels the internal en- 
ergy is 

U = (H(T-r)) Trt 

= J2 I drV(T\r,p)V(r\t)V(t)H(T;r), (26) 

*r 4- •* 



where 



V(t\t,P) oc e 



-0H 



(27) 



is the thermal Gibbs probability at inverse temperature (3 
which obeys V(r\r, = 1) = V{t\t). Since under such a 
gauge transformation the Hamiltonian remains invariant, 



(21) we have H(t;r) = W(l;tr), where tr = (t 1 



and 1 is an M-dimensional vector with all entries equal 
to 1. Therefore, one can write 



[/ = W dr 



V(r\r,P)V(T[ 
V(t\P) 



-V(r\t)V(t)H(l;TT) 



(28) 

Gauging the variables rr — > r, reorganising the terms 
and taking (3 — 1, we finally get 



U = J drP(r|l)W(l;r). 



(29) 



The meaning of this is that, for a gauge invariant chan- 
nel of the type described above (which includes the sym- 
metric channels), the internal energy is independent of 
the configuration. In special cases, as can be found in [2], 
the gauge symmetry allows an analytical expression to be 
found. The same method can be used to prove that the 
probability distribution for the magnetisation is equal to 
the probability distribution for the two-point correlations 
in Nishimori's temperature, which indicates the absence 
of a spin glass phase and no replica symmetry breaking. 
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FIG. 2: Schematic drawing of the relay-without-delay (RWD) 
setup to be analysed. 



|Sh, union's Limit 



- Relay-without-Delay 
■ Gaussian Channel 



A. The Relay Channel 

In order to compare our results with those of [l6|, we 
analyse the RWD for the setup sketched in Fig. [2] with 
a\ = -qa 2 , a = h = 1 and c% = (1 + a 2 )- 1 / 2 . The 
potential is then given by 



*(r,r) = 




+ e -{r-r+l?/2^ erfc ( + \ i (30) 

where erfc (a;) is the complementary error function 

erfc (re) = —= \ e~ y2 dy. (31) 

V 7 *" Jx 

For these values of noise and gains, the capacity of this 
channel as derived in [l6[ is 

C = \ log 2 (l + . (32) 

The numerical results for the overlap between the re- 
trieved and the original codewords, obtained by solving 
recursively equations (fT7|l . are given in Fig. [3] for K = 4, 
C = 3, j3 = 1 and r\ = 0.1. Shannon's limit is indicated 
by the vertical dashed line and corresponds to a noise 
level a 2 « 8.79. The dashed curve shows the overlap 
for a simple Gaussian channel with noise level a 2 and 
the continuous one shows the overlap for the RWD. The 
improvement in the practical limit for error-free commu- 
nication, indicated by the highest noise level for which 
d = 1 is clear. However, the distance between the dynam- 
ical transition threshold a 2 , 2.22 and Shannon's limit 
for the channel is greater than in the case of the simple 
Gaussian channel (for numerical results for the Gaussian 
channel see [20j). Numerical calculations point to the 
expected result that decreasing the noise level from the 
source to the relay brings a 2 , closer to Shannon's limit. 
However, one must remember that the relay strategy ex- 
amined does not use the full potential of the relay and 
the additional information embedded in the LDPC codes. 
We expect that a LDPC decoding in the relay will im- 
prove the communication performance and currently fo- 
cus on the analysis of this scenario. 



FIG. 3: The overlap of retrieved and original codewords for 
the RWD channel with a\ = 0.1a 2 is given by the continu- 
ous curve. The dashed curve shows the same for the simple 
Gaussian channel. The vertical line indicate Shannon's limit 
for the RWD as calculated by El-Gamal and Hassanpour [la ]. 



We can also see in the plot that, as the noise level 
increases, the channel becomes closer to the Gaussian 
channel. This is just a consequence of the fact that, for 
high noise level, the additional information provided by 
the relay becomes negligible as both users decode the 
message poorly. 

Figure [?] shows the entropy and the free and internal 
energies for the same values as in Fig. [3J At the dynam- 
ical transition point, where practical perfect decoding 
becomes unfeasible, the entropy becomes negative, indi- 
cating the emergence of subdominant metastable states 
that can be further explored using the replica symmetry 
breaking ansatz. Between this point and the thermody- 
namical transition point, where the entropy is positive 
again, the dominant state is still ferromagnetic but the 
population dynamics algorithm used to solve the sad- 
dle point equations becomes trapped in a local minimum 
with free-energy higher than the ferromagnetic one. Due 
to the equality between the internal energy and the fer- 
romagnetic free-energy, the point where the entropy be- 
comes positive again is also the point where both energies 
cross in the bottom graph. 

In Fig. [5] we plot the dynamical and thermodynamical 
transition noise levels against 77, the ratio between the 
relay and final receiver noise levels. We see that the dy- 
namical and thermodynamical transition points decrease 
with rj but became closer to each other, stabilising at 
asymptotic values that match those of the simple AWGN 
channel values as the relay contribution becomes mean- 
ingless. 

Although the capacity for the RWD is known only in 
special cases, its upper bound can be higher than in the 
case of the CRC. In order to verify it for the LDPC-based 
framework we analyse in this paper, we now use a setup 
equivalent to the one studied in [151 ] and shown in Fig. [S] 
where a = c\ = b± = 1 and a 2 = Act 2 . The capacity of 
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FIG. 4: Entropy and energies for the RWD channel at Nishi- 
mori's temperature. The upper graph shows the entropy 
which is given by the difference between the internal energy 
and the free energy depicted in the bottom graph. 
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FIG. 5: The continuous/dashed line shows the dynam- 
ical/thermodynamical transition noise levels of the RWD 
against r), the ratio between relay and final receiver noise lev- 
els. The upper horizontal line corresponds to Shannon's limit 
a 2 « 8.79. 



FIG. 7: The continuous/dashed line shows the dynam- 
ical/thermodynamical transition noise levels of the RWD 
against A in the setup of Fig. [6] The continuous line without 
marked symbols is Shannon's limit for a CRC with the same 
noise levels and transmission powers. 

the CRC in this case is 

f I l0 S2 (1 + 7^). *>1, 
C={ (33) 

[ 3 lo g2(l + TT+wj' A<L 

In Fig. [7J we compare the dynamical and thermody- 
namical threshold noise levels of a RWD with Shannon's 
limit for the CRC, both with the setup described above, 
for different values of A, the ratio between the noise levels 
applied at the transmission and reception points. 

We can see that, although the practical decoding line 
(dynamical transition) falls below Shannon's limit for all 
calculated values, the thermodynamical transition goes 
above it for the CRC case at higher values of A. Figure [7] 
shows that the capacity of the RWD is indeed higher 
than the CRC for the case studied and quantifies the 
gain in allowing the message sent by the relay to depend 
on the current transmitted symbol (which is excluded in 
the CRC). Although allowing this instantaneous depen- 
dence would at first sight seem just a small modification 
which is insignificant in the infinite block length limit, 
it indeed gives relevant extra information which facili- 
tates more efficient retrieval at the final receiver. The 
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insight gained is that for the RWD and large A, the relay 
transmission t± is correlated with the original codeword 
t M , which is not the case in the CRC allows; this allows 
for an improvement in the information extraction at the 
receiver. 



B. Large Relay Array 

Now, we will use the central limit theorem to obtain the 
result for large L in the relay array setup given by Fig. [I] 
As the relay messages are correlated and to guarantee 
that the quantities have the same order, we introduce a 
1/L scaling in the summation over relay messages. The 
potential for this model becomes 



*(t, r) = J < dri exp 



x exp 



^W( r " r ^? sgnri ) 



(34) 



where we assumed, for simplicity, a = hi = Cj = 1. 

For L ^> 1, the central limit theorem amounts to a 
modification in the distribution of the variable r given 
by 



where 



Y[dnV(j 




V27 



: exp 



exp 



(r- 1 -uof 

{n - if ' 

2K? + <r 2 ) 



(35) 

(36) 

(37) 
(38) 



a-n Thermodynamical Transition 
O — o Dynamical Transition 




L (Relay units) 



FIG. 8: Dynamical and thermodynamical transition points 
for many relays. The exact formula is used to calculate the 
points L = 1, 2, 3, 4, 5. The horizontal lines represent the large 
L limit. 



the limit L — > oo, this distribution becomes a delta func- 
tion centred at the error function and therefore 

V(r)=Af(l + ed(l/V2aXafj. (42) 
Accordingly, the potential becomes 



1 



*( r) r)=expj-^ 



t — erf \t I \/2a, 



(43) 



Figure [8] compares the dynamical and thermodynami- 
cal transition points for L = 1, 2, 3, 4, 5 calculated by the 
exact formula and the result obtained by the approxima- 
tion for large L. Again, we consider the case of K = 4, 
C = 3 and 13 = 1. 

It is clear from Fig. [5] that already at L = 5, both 
dynamical and thermodynamical transition points ap- 
proach the large L limit solution, thus making this ap- 
proximation attractive already for low L values. 



with 



\ i=l i=l / 

(39) 

For simplicity, we consider the case where the noise 
level is the same for all relays of = o\ and define 



^ 2 — ~ 2 
o„ = o. 



2 



~2 — „2 
Of = O 



'ituqi u f = u TW o- (40) 
Then the corresponding distribution for to 



V(uj) =Af( erf (l/%/2o r ),ierfc 2 (l/V2o, 



(41) 



Consequently the contribution for the final noise level 
coming from the relay transmission decreases as L" 1 . In 



V. CONCLUSIONS 

In this work we analysed the behaviour of relay arrays 
using methods of statistical mechanics. These networks 
are of growing significance due to the increase of multi- 
user, mobile and distributed communication systems. 

We found an analytical solution for the relay-without- 
delay (RWD) channel given by the RS ansatz, which due 
to the gauge symmetry of the channel, is exact at Nishi- 
mori's temperature that correspond to a choice of the cor- 
rect prior within the Bayesian framework. We showed the 
level of improvement with respect to a simple Gaussian 
channel without relaying which, even for the naive re- 
lay strategy of Demodulate-and-Forward analysed here, 
is significant. 

We compared the RWD dynamical and thermodynam- 
ical transition points for different noise ratios between the 



8 



relay and the direct channel; and found that although 
these points become far from Shannon's limit, the differ- 
ence between the dynamical and the thermodynamical 
transition decreases. The relevance of the relay is clearly 
decreasing as its noise level increases as the level of ad- 
ditional information it conveys diminishes. 

We also were able to compare the RWD case to the 
classical relay channel (CRC) for different noise ratios 
between the relay and the direct channel. We found that 
the capacity of the RWD is higher than the CRC for a 
high relay noise, showing the significance of the extra 
information conveyed by the relay on the current trans- 
mitted symbol, which is absent in the CRC framework. 

The performance of a large array of relays was anal- 
ysed and compared against results obtained for a small 
number of units. The result obtained are consistent and 
indicate that this useful approximation provides accurate 
results already for a small number of units. For a large 



array, we also found that the increase in noise tolerance 
levels off. 

We have demonstrated the usefulness of methods 
adopted from statistical physics for analysing multi-user 
communication systems. While we have concentrated 
on limited scenarios of relay channels, we believe that 
these methods hold a promising alternative to the infor- 
mation theory methodology which, in general, has not 
been successful in dealing with multi-user communica- 
tion systems. The study of different relay channels and 
other multi-user communication networks is underway. 
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